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I. INTRODUCTION 



The strong influence of the electron dynamics provides the possibility of controlling the expansion 
of laser-produced plasmas by appropriately shaping the laser pulse. A simple irradiation scheme is 
proposed to tailor the explosion of large deuterium clusters, inducing the formation of shock struc- 
tures, capable of driving nuclear fusion reactions. Such a scenario has been thoroughly investigated, 
resorting to two- and three-dimensional particle-in-cell simulations. Furthermore, the intricate dy- 
namics of ions and electrons during the collisionless expansion of spherical nanoplasmas has been 
1 analyzed in detail using a self-consistent ergodic-kinetic model. This study clarifies the transition 

O |' from hydrodynamic-like to Coulomb-explosion regimes. 

PACS numbers: 36.40.Gk, 52.38. Kd, 52.65-y 

11: 

The rapid developments in laser technology in recent years, both in the IR [1| and VUV/X frequency range 0], 
c/3 , disclosed new research realms in the field of radiation-matter interactions. Among these, the interaction of ultraintense 
lasers with jets of atomic or molecular clusters has become a central research topic Experiments with laser-cluster 
*-£h I interactions have demonstrated tabletop nuclear fusion 0, opening the way to compact neutron sources; recent 
studies also suggest the use of the laser-cluster scheme for laboratory investigations of nucleosynthesis reactions, 
relevant in astrophysical studies Q . 

A clustered gas beam can be obtained through adiabatic expansion of a dense gas jet into a vacuum 0, H], with 
the gas backing pressure and the temperature determining the mean cluster size (typically, in the range 1 nm - 1 
/im). Such a medium can be regarded as a sparse distribution of tiny solid- like targets, a peculiar configuration that 
allows for a deep penetration of the laser radiation and a strong laser-matter coupling with many individual, overdense 
targets: this guarantees an extremely efficient absorption of the radiation (nearly the entire laser energy deposited 
f*^ . within a few millimeters propagation length 0). 

f* 1 **. ' When an ultraintense laser beam hits a cluster, the leading edge of the pulse promptly ionizes the neutral atoms 
O ! (cf. Ref. [3 for a detailed analysis of the concurring ionization mechanism in different laser/cluster configurations), 
*"^5 ' forming a dense distribution of ionized matter (usually called a nanoplasma) 0, [ll[ . Then, the electrons are further 
heated by the main part of the pulse and the plasma starts to expand, leading to efficient ion acceleration, as first 
r>3 ' predicted by Dawson [J. M. Dawson, "On the production of plasmas by giant laser pulses", Phys. Fluids 7 pp. 981- 
£>"\ 987 (1964)]. For small, low-Z clusters (e.g., nm-sized Hydrogen or Deuterium clusters) exposed to extremely intense 
laser radiation, the electrons can even be completely stripped from the cluster in a few optical cycles, leading to the 
Coulomb explosion (CE) of the remaining bare-ion distribution. With IR lasers, this happens when the excursion 
length of the electrons \c/u>i arcsin (ao/(l + a 2 ))), with u>i and eto as central frequency and normalized vector potential 
of the laser, respectively] is much higher than both the electron skin depth (c/w pe , to pe being the electron plasma 
frequency) and the initial radius of the cluster, i?o fl2| . In more general situations, when a relevant fraction of the 
electrons is bound to the cluster, the expansion process is strongly dependent on the self-consistent dynamics of ions 
and trapped electrons: when increasing the laser energy or when lowering the cluster size, the expansion conditions 
vary smoothly from quasi-neutral, hydrodynamic-like regimes to pure CE regimes, as revealed b y p article-in-cell (PIC) 
simulations [l2|, HH of laser-cluster interactions and detailed kinetic analysis of the expansion [1J] . This suggests the 
possibility of exploiting the electron dynamics in order to achieve control over the expansion, which can be done by 
appropriately shaping the laser pulses. For example, large-scale shock shells (multi-branched structures in the ion 
phase space, previously predicted by Kaplan, Dubetsky, and Shkolnikov [l5| for nonuniform CEs) can be driven and 
controlled by using a proper sequence of two laser pulses having different intensitie s 11611: t his opens the way toward 
intracluster fusion reactions within large Deuterium or Deuterium- Tritium clusters [la , Il7j | . 

The present paper is organized as follows. After briefly introducing the concept of shock shells in a Coulomb 
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explosion (Sec. Ill A[) . the description of a double-pump irradiation scheme (Sec. Ill Bp for inducing and controlling 
shock shells capable of driving intracluster nuclear reactions (Sec. Ill CP is given, whereas Sec. Ill Dl investigates the 
possibility of controlling expansions also with VUV/X-ray beams. Section Mil is devoted to the kinetic analysis of 
spherical plasma expansions, performed with a novel ergodic- kinetic model (Sec. IIII A|) : selected results obtained with 
this model are presented in Sec. IIII Bl 



II. CONTROLLED SHOCKS SHELLS 



A. Shock shells in Coulomb explosions 



The explosion dynamics of a spherical distribution of cold ions is described by the equation for the radial motion 
of the ions. 
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where ri(ro, t) denotes the trajectory of an ion starting at t = from radius ro with zero radial velocity, M is the ion 
mass, qi = Ze is the ion charge, and N{(r, t) is the number of ions enveloped at time t by a sphere of radius r. When 
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FIG. 1: Initial ion density (blue) and radial electric field (light red) for the case of (a) a uniform, step-like density and of (b) 
a nonuniform, smoothly-decreasing density, defined as nio/ [l + (ro/Ro) 6 ] 3 ^ 2 , Ro being the radius of the equivalent uniform 
sphere. 



starting from a uniform ion density (equal to n; for r ^ R ), the repulsive electric field grows linearly for r < R Q , 
reaching its maximum at the outer boundary (Fig. [TJt), and causing the outer ions to be always faster than the inner 
ones. In such a particular situation, ions never overtake each other, iVj(r, t) is conserved along the ion trajectories (so 
is the total energy of each ion), and Eq. (fTJ) can be integrated analytically [l5|, [l^: assuming all ions to be initially 
at rest, the expansion velocity Vi(ro,t) is 
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where uj p i = y/4irqfn iQ /M is the (initial) ion plasma frequency, whereas its radial trajectory is given by 

VI - 1) + log ( v 7 ^ + V^i) = yfwpit, (3) 

where £j = r;/ro is the expansion factor. Since the right-hand side of Eq. ([3]) does not depend on ro, £i is the 
same for all ions, being a function of time only, i.e., £ = As a consequence, the i> — r phase-space profile is 

always a straight line [having equation v = v/2/3(£ — l)/£ 3 Wpir] and the ion density stays always uniform (evolving 

as n £~ 3 ). Correspondingly, the ion energy distribution is 3/2(^0/9?) \A/ £ ce(1 — l/£)~ 3 / 2 , where 6ce = qfNio/Ro 
is the asymptotic cutoff energy. Such solutions could be obtained equivalently from cold-fluid equations (v being 
now the fluid velocity). However, when nonuniform density profiles are considered, the kinetic aspects of the process 
immediately reveal themselves and the features of the expansion change dramatically, as first predicted in Ref . [l5| , 
and explored in self-consistent simulations in Refs. (ltL fl7j. In fact, if the initial density is a decreasing function of 
r, the repulsive Coulomb field reaches its maximum within the ion sphere [HI, [l5| (cf. Fig. Wp), causing some inner 
ions to accelerate more than the outer ions in front of them: this ultimately leads to overtaking between ions and the 
formation of a shock shell. In this situation, the ion trajectories arc no longer independent of one another, N{(r,t) 
is not conserved along the ion trajectories (nor is the total energy of each ion), and Eq. JT]) no longer admits the 
analytical solutions reported above. Figure [5] shows a typical nonuniform expansion (similar to that illustrated in Ref. 
[lH . Fig. 1), starting from the smoothly-decreasing profile of Fig. [TJd: the phase-space profile starts bending on the 
right until it becomes multivalued, and a pronounced three-branched shock shell forms (cf. Fig. [5^). As predicted 
by the theory [l^|, vertical-tangent points in the phase-space profile correspond to a singularity in the ion density 
(cf. Fig. (2b) (such singularity is sometimes called a caustic using a term borrowed from astrophysics (2(|), while 
horizontal-tangent points correspond to a singularity in the energy distribution (cf. Fig. [2}:); this appears at the very 
beginning of the expansion and it is sometimes called the shock predictor [l5| . since it indicates the presence of a 
maximum in the velocity profile. As the expansion continues, the shock shell widens radially, involving an ever larger 
portion of the ion cloud. The physical interest of such pronounced, large-scale shock shells resides in the appearance 
of large relative velocities within a single exploding cluster, which can lead to energetic ion-ion collisions, thus opening 
the way to intracluster nuclear fusion reactions [H, [IB, E3 ( c -g-; m the case of D or D-T clusters). 
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B. The double-pump technique 



An efficient strategy for producing large-scale shock shells in a controlled fashion is the combination of different 
expansion regimes. In fact, if a slow, hydrodynamic-like expansion is exploited to provide a decreasing density profile, 
then a subsequent, abrupt CE naturally leads to the formation of a shock shell, according to the dynamics described 
in Sec. Ill Al Such a slow-expansion/abrupt-explosion dynamics can be obtained rather simply, by irradiating a cluster 
with two sequential pulses having different intensities, namely, a weaker pulse followed (after a proper time delay) 
by a much stronger one, as sketched in Fig. [3l The first pulse must be intense enough to ionize the atoms, creating 
a nanoplasma, but no so intense as to expel a significant fraction of the electrons from the cluster core: in these 
conditions, a slow expansion takes place, in a quasineutral, hydrodynamic-like regime, and a smoothly decreasing 
density profile is formed. The second pulse must be extremely intense, in order to strip all the electrons from the 
cluster, thus driving a CE. The effectiveness of this double-pump technique has been demonstrated [HI, E3] by 




FIG. 3: Sketch of the double-pump scheme: the target cluster is first irradiated with a low-intensity laser and, after a delay 
At, with a high-intensity laser. 




FIG. 4: Ion density at (a) t = 170 fs, (b) t = 187 fs, (c) t = 206 fs, and (d) t = 237 fs. A lfim. x l^tm computational domain 
has been used; a 840 x 840 uniform spatial grid, and 1.2 x 10 6 particles per species have been employed. 

resorting to 2D and 3D PIC simulations performed using the OSIRIS 2.0 framework [2l[, closely matching realistic 
physical scenarios. The results presented in Figs. 2] and [5] refer to a 2D simulation, in which a deuterium cluster 
(initial radius Rq = 32 nm, density no = 4.56 x 10 22 cm -3 ) is first irradiated by a relatively weak laser pulse (peak 
intensity 1 x 10 16 W/cm 2 , duration 35 fs, wavelength 820 nm) and, after a 170 fs delay, it is irradiated by a second, 
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FIG. 5: Ion phase-space profile at times tO = 170 fs, tl = 187 fs, t2 — 206 fs, t3 = 225 fs, and t4 = 237 fs, for particles contained 
in an angle A9 — 0.1 rad around the propagation direction x (blue markers) and the polarization direction y (light-red markers). 

extremely intense pulse (peak intensity 2.5 x 10 19 W/cm 2 , duration 20 fs, wavelength 820 nm); both pulses propagate 
in the x direction, being linearly polarized along the y direction. Figures [5] and [5] show the evolution of the ion density 
distribution and the ion phase-space profile, respectively, during the CE stage. When the second pulse hits the plasma 
sphere (cf. Fig. S^,), the outer ions are expanding slowly, under the pressure of the hot electrons, while the cluster core 
is still dense, with the inner ions still at rest (cf. Fig. [5]). As a result, when the electrons are suddenly swept away by 
the second laser, the inner ions experience a much higher repulsive force than the outer ions do, so that they rapidly 
overrun them, leading to the formation of a large, pronounced shock shell, clearly visible in both ion density (Fig. 
HJi) and phase space (Fig. [5]). In this example, the relative velocities between ions belonging to different branches in 
the shock shell would allow for collision energies as high as 200 keV, which is enough for efficient DD fusion to occur. 

C. Fusion reactions within single exploding clusters 

The ability of generating pronounced shock shells involving high relative velocities within single exploding clusters 
makes the phenomenon attractive as a possible way to induce nuclear reactions. As already pointed out in Ref. 
[lil ]. the rates for intracluster reactions can be much higher than those for intercluster reactions, because the typical 
densities within a shock shell (< nio ~ 10 22 cm -3 ) can be much higher than those within the hot plasma filament 
resulting from the exploded clusters, whose density is ~ 10 19 cm -3 . However, it is crucial to consider also that a 
shock shell stays appreciably dense only for a very brief time (~ 10 fs) , much shorter than the typical disassembly 
time of the plasma filament (~ 10 — 100 ps). In fact, as shown in Ref. [ijj, the intracluster-reaction rates exhibit 
a sharp, time-resolved peak right after the shock-shell formation. In order to compute the actual number of shock- 
driven, intracluster fusion reactions, it is necessary to sum up all possible contributions from collisions between ions 
belonging to different velocity branches. The number of reactions per unit time and volume, 1Z, is given by 

7i = y^ n h (r)n k (r)a (\v h - v k \) \v h - v k \, (4) 

h<k 

where a is the cross section for DD fusion, while and Vh indicate, respectively, the ion density and velocity on 
the hth branch. The number of intracluster fusion reactions, Af, is obtained by integrating 1Z over time and space, 
as Af = 47T 7^°° / r Rsh lZr 2 drdt, where r s h and i? s h are the shock shell boundaries, whereas i s h is the shock formation 
time. Clearly, Af is strongly dependent on the dynamics of the shock shell, which can be tailored by varying the 
delay between the two laser pulses as well as their intensities. A detailed analysis of the influence of the double-pump 
parameters on the intracluster fusion rate (cf. Ref. [l7| ) reveals that optimal delay/intensity combinations exist that 
maximize Af. For very large clusters (i?o ^ 100 nm) and optimal double-pump configuration, the intracluster reaction 
yield can become comparable with the interluster neutron yield, with ~ 10% of the fusion reactions arising from 
intracluster collisions 1 171 . Since intercluster reactions occur on a much longer time scale than intracluster ones (1Z is 
sharply peaked in time [17j]). a well tuned double-pump experiment should provide a clear signature for the occurrence 
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of intracluster, shock-driven fusion reactions, in the form of a timc-resolvcd burst of fusion neutrons anticipating the 
usual bulk of fusion neutrons due to intcrcluster reactions. 



D. Expansion control with intense VUV/X-ray beams 

The possible use of the double-pump technique also with VUV/X-ray radiation is currently being investigated. With 
very short wavelengths, the expulsion of electrons from the cluster, which is necessary for a Coulomb explosion to 
occur, depends essentially on the ionization dynamics and much less on the electron motion in the radiation field [since 
now the electron excursion length (~ c/uji) is very small compared to the typical cluster radius]. Thus, the full details 
of the ionization processes involved in this situation (e.g., photoionization, Auger decay, collisional ionization, and 
strong field ionization) must be considered. For example, when the main ionization process is the photoelectric effect 
(as in D clusters), the mean energy acquired by the electrons is approximately equal to the photon energy, e 7 = hv. A 
Coulomb explosion can occur if the number of photons hitting the cluster is sufficient to ionize most of the atoms, and 
e 7 is higher than the electrostatic potential well generated by the ions. This suggests that a double-pump technique 
can still be employed to drive a shock shell, even though the physical mechanisms it relies on are different. In fact, 
the goal of the first, weak pulse is now to ionize a small fraction of the atoms in the cluster, in order to drive a slow 
expansion, whereas the second, strong pulse must ionize all the atoms remaining in the cluster core, causing a sudden 
explosion. In these conditions, a shock shell forms because some of the ions created by the second pulse overtake 
the slowly expanding ions, provided that appropriate radiation intensities and delay between pulses are employed. In 




FIG. 6: Ion phase-space profile at times t\ — 160 fs, t-2 = 250 fs, tz — 390 fs, and ti = 460 fs. The computational domain is 
a cube with a 1.2 fj,m side (256x256x256 cells), and 125 particles per cell per species have been used (in the region initially 
occupied by the cluster). 

order to test the effectiveness of this technique, PIC simulations have been used [2l[ , which describe the evolution 
of the cluster given the initial distributions for the electron and ion populations. This approach is valid whenever 
the ionization time is short with respect to the expansion time scale, which is usually a good approximation for the 
intense light sources available in the VUV/X-ray energy range. The formation of a pronounced shock shell, obtained 
using this technique, is shown in Fig. [5J referring to an ideal situation in which it is assumed that a deuterium cluster 
(initial radius Rq = 30 nm, density no = 10 22 cm~ 3 ) undergoes ionization processes which lead to the formation of a 
population of Ao/3 electrons (where Nq is the number of initially neutral atoms) with 12.4 keV energy, followed by a 
second population of 2Nq/3 electrons with equal energy, after a 150 fs delay. We observe that at ultra high intensities 
in XUV/X-ray sources the ionization dynamics is poorly understood since strong field ionization and stabilization 
against this mechanism may come into play (22j. The balance between the photoelectric, strong field ionization and 
stabilization processes must be analyzed in order to obtain the correct picture for the energy distribution of the 
electrons (which actually determines the dynamics of the explosion). The analysis of double-pump scenarios, taking 
into account the full dynamics of ionization, with beam configurations available at LCLS and XFEL and with different 
cluster constituents will be presented in future publications. 
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III. KINETIC ANALYSIS OF SPHERICAL PLASMA EXPANSIONS 



The possibility of controlling the expansion by acting on the amount of energy transferred from the radiation to 
the electrons stresses the need for a deep understanding of the physics of the electron-driven collisionlcss expansion of 
spherical, nm-sized plasmas, composed of cold ions and hot electrons. This is also crucial for particular applications, 
such as the biomolecular imaging with ultrashort X-ray pulses [23j . where sample damage before the imaging time 
must be avoided. 

Since the first experiments on laser-cluster interactions were performed || , it has been debated [ill . [2H whether 
the cluster expansions were driven mainly by the pressure of the hot electrons (hydrodynamic-like regime) or by the 
strong repulsive field which forms when a large fraction of the electrons leave the ion core (CE). Actually, these two 
expansion regimes are just two opposite limits resulting from the same physical process, namely the formation of a 
repulsive electric field due to energetic electrons trying to escape from the cold-ion core. When the kinetic energy of 
the electrons is very small compared to the electrostatic energy stored in the ion core, the electric field is localized at 
the ion front [lH l25l. |26| . the whole cluster being approximately neutral, and the expansion is hydrodynamic-like; on 
the contrary, when the kinetic energy of the electrons is much higher than the ion electrostatic energy, most of the 
electrons are free to escape and the CE of the ion sphere takes place. However, in a general scenario, the dynamics 
will be a mixture of the phenomenology of the two limits, and a kinetic description is necessary in order to capture 
the detailed features of each particular regime. 



A. Ergodic model 

In order to analyze accurately the dynamics of the expansion in any intermediate situation, we have developed 
a new kinetic model, which allows for a clear identification of the transition from the hydrodynamic expansion to 
the CE regime, thus setting the range of validity for the CE approximation [l4j|. Based on the assumptions that 
the electrons are nonrelativistic and that the time scale for their motion is much shorter than the time scale for the 
ion motion, the model describes the expansion dynamics self-consistently by following the radial motion of the ions, 
whereas treating the evolution of the electron distribution as a sequence of equilibrium configurations, represented 
by stationary solutions of the Vlasov equation for the electrons. In the case of ideal, spherical symmetry, the general 
stationary solution is a function of the Hamiltonian and of the absolute value of the angular momentum, £, which are 
constants of motion of the system. However, and by noticing that possible perturbations to the spherical symmetry 
would rapidly break the ^-conservation, while leaving the ^-conservation unaltered, the distribution function in 
equilibrium can be assumed to depend on Jtf only, i.e., the system can be assumed to be ergodic. Accordingly, the 
electron density in equilibrium with the electrostatic potential $ can be written as 

ne(r) = i ^|p(e)^(r,e;{$})de J (5) 

where p is the energy distribution of the electrons and 

&(r,e; {*}) = - ( ' ± (6) 

r' [e + e$ (/)] 2 dr' 

is the probability of finding an electron having total energy e at the radial position r. The self-consistent potential <f> 
satisfies the Poisson equation V 2 $ = 4ire(n c — rij) , which is nonlinear because n depends on $. In order to obtain a 
closed set of equations describing the plasma expansion, one still needs a rule to update the equilibrium configuration 
when the ion motion is considered. Since the variations of $ induced by the ion motion are slow when compared with 
the time scale for the electron motion, the theory of adiabatic invariants for time varying Hamiltonians is used. Here, 
the appropriate invariant is the ergodic invariant (cf., for example, Rcf. [27) |) 

J ( e ) = ^^ TO 3/2 j (£ + e$) l r 2 dr ^ (7) 

which is the phase-space volume enclosed by the surface |mv 2 — e<!> = e. The conservation of X means that the energy 
variation can be expressed as the ensemble average of the potential-energy variation, thus providing the desired rule 
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for the evolution of the electron energy. Finally, a closed set of equations can be written in the form 




47re (n e — Zn\) 



(8d) 



(8b) 



(8a) 



(8c) 



(8e) 



where Eqs. (|8a|) and ([5c| . determining the evolution of the ion trajectories n(ro, t) and the electron energies e(eo, t) (eq 
being the initial energy), are coupled with the nonlinear Poisson equation, which provides the self-consistent <I>(r, t) 
(here, $ is set to zero at infinity, so that e < for trapped electrons). In Eq. (|8dj) . the relation p(e,t)de = po(eo)deo 
has been used, po(eo) being the initial energy distribution of the electrons. Equation (|8cj) expresses the ion density 
in terms of the Lagrangian coordinate ro, under the hypothesis of no ion overtaking (cVj/cVo ^ 0) [l5| (when shock 
shells are present, a more complicated expression has to be used (HI)- The solution of Eqs. (8) requires the knowledge 
of riio(ro) and po(eo)- Depending on the physical problem considered, the determination of pa(eo) can be nontrivial. 
This is precisely the case when considering expansions starting from an initially neutral distribution (n c = Zn\, and 
$ = everywhere) with hot Maxwellian electrons, having temperature To, because such configuration is far from 
equilibrium. In that case, the plasma expansion can be regarded as a sequence of two distinct processes: first, a 
rapid expansion of the electrons is observed, which leads to a VP equilibrium configuration before the ions move 
appreciably; second, a slow expansion of the plasma bulk is observed. While the latter stage can be certainly analyzed 
using Eqs. (8), a different procedure must be used to determine the equilibrium configuration following the first stage 
[in order to obtain the correct form of po(eo)]- To this purpose, the actual charging transient described by the full 
VP model is replaced by a virtual charging transient, in which an external potential barrier, initially confining the 
electrons, is gradually moved from Rq to infinity, with a scries of small radial displacements. Each time the barrier 
is moved farther (from R w to R w + SR W ), the new self-consistent potential ($ + S$>) and electron energy (e + Se) 
have to be calculated. In order to mimic an expansion into a vacuum (which the real transient is), direct energy 
exchanges between the electrons and the expanding barrier (expansion work) must be avoided. This is accomplished 
by expressing the energy variation Se as 



Although formally identical to Eq. (|5e|) . Eq. © is not equivalent to the conservation of the ergodic invariant (which 
would bring on an extra term in the expression for Se, accounting for the expansion work). The variations 5$ and 
Se are determined by simultaneously solving Eqs. (|8b[) and ([9]) with a suitable iterative scheme. Once coupled with 
this procedure, which provides the initial conditions po(eo), Eqs. (8) constitute a self-consistent model describing 
the collisionlcss expansion of spherical nanoplasmas, starting from an initially neutral distribution composed of cold 
ions and hot, Maxwellian electrons. This model has proved to be suitable for a detailed analysis of the expansion 
mechanism in a wide range of conditions, since the dynamics can be determined by the single dimcnsionlcss parameter 
To = ZfeeTo/ecE = 3A^/i?o (^d being the initial Debyc length for the electrons), which accounts for both the initial 
electron temperature and the cluster parameters (size and density). 



In the following, a detailed analysis of the expansion dynamics is presented for two reference cases, (a) T = 7.2 x 10~ 3 
and (b) To = 7.2 x 10~ 2 (for the cluster size and the density employed in Sec. Ill Bl these values correspond to T) = 2 
keV and 20 keV, respectively), also illustrating the whole-range To-dcpendcnce of the key expansion parameters. 

Figures [7] and [5] show the self-consistent equilibrium configuration of the electrons after the initial charging transient 
(comparisons with reference results obtained with OSIRIS are also shown, which confirm the validity of the barrier 
method, outlined before, to obtain the conditions after the charging transient). In Fig. [7| the electron density is 
plotted, along with the corresponding electric field: the positive charge buildup at the ion front, AQ, is in (a) 12.5% 




(9) 



B. 



Selected results 
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FIG. 7: Equilibrium electron density (blue) and electric field (light red), after the charging transient with immobile ions, for 
(a) Tb = 7.2 x 10 -3 and (b) To = 7.2 x 10~ 2 . Lines refer to the ergodic model and markers to PIC simulations. 




FIG. 8: Equilibrium electron energy distribution po(eo), after the charging transient with immobile ions, for (a) To = 7.2 x 10 
and (b) To = 7.2 x 10 -2 . Lines refer to the ergodic model and markers to OSIRIS simulations. 



and in (b) 38% of the total ionic charge eNo. Figure [8] shows the equilibrium energy distribution po(eo), to be used 
as initial condition for the bulk expansion. The To-depcndence of AQ and of the mean kinetic energy of the trapped 
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FIG. 9: Charge buildup AQ (blue) and mean kinetic energy of trapped electrons, £ , (light red) as functions of To, after the 
charging transient with immobile ions. Circles refer to the ergodic model, solid lines to the fit laws of Eqs. (|10p and (|ll|l . 
Dashed lines show the power-law behavior of AQ and £ for To -C 1. 

electrons, £, is displayed in Fig. [9] Simple, accurate fits exist for these quantities, as 




37^ = 1 -^"3.35 (l.86f 1/2 ), (11) 
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where .7>(a0 = x/(l + a;' 1 ) 17 ''- For T < 1 (i.e. A < i?o), Eq. JTU]) reduces to AQ/Q ~ v/o/eTo , thus recovering 
the theoretical results for planar expansions [26|, [28|, while Eq. (fTTj) reduces to £/(§ /cb^o) ~ 1 — 1.86 . The 




Time [oi ,1 Time [u 1 ] 

L pi J L pi 1 

FIG. 10: Evolution of the ion-front trajectory (panels a and b), of the electronic charge enveloped by the ion front, eiVo — AQ, 
(panels c and d, blue) and of the mean kinetic energy of trapped electrons, £, (panels c and d, light red), for (a,c) To = 7.2 x lCP' 5 
and (b,d) T = 7.2 x 10" 2 . 

time evolution of AQ and £ during the bulk expansion is illustrated in Fig. 1101 as the ions expand, gaining kinetic 
energy, the electrons cool down and the charge buildup decreases (until a ballistic regime is reached for both species 
[29l]). Such behavior of the electrons strongly affects the ion dynamics and their resulting energy spectrum, since the 
asymptotic energy f m of an ion is given by 

£qoM = g(r ,0) + r°° 1 dq (n(ro, t),t) 
Ze r J Q n(r ,t) dt 

where q(r,t) is the net charge buildup enveloped by a sphere of radius r at time t. The integral term (vanishing 
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FIG. 11: Asymptotic ion energy spectra for different values of To, compared with the theoretical asymptotic spectrum for the 
CE case. 

for a CE) accounts for the energy loss due to the decrease of the positive charge buildup experienced by each ion 
along its trajectory. Figure [TT] shows the asymptotic ion spectrum for different values of To. For T < 0.5, profound 
differences exist with respect to the CE case, with the distribution exhibiting a local maximum far from the cutoff 
energy. The spectrum is monotonic only for To > 0.5 (a condition already close to the CE regime, with cutoff energy 
above 0.7£ce)- Therefore, the transition value To = 0.5 can be considered as the lower bound for the validity of the 
CE approximation. The cutoff energy e max and the location of the maximum, e pca k, are plotted in Fig. [T^jas functions 
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of T . The cutoff energy can be accurately fitted as 

£max = ^1.43 (2.28 £ C E, (13) 

which, for To <C 1, reduces to e max ~ 2.28Tq^ 4 £ce, while, for To < 0.5, e poa k exhibits the simple power-law behavior 
fpeak = 0.3TQ' 9 ecE- These formulae provide accurate scaling laws, valid for any combination of i?o, wp, and To (as 
long as relativistic effects can be neglected) , which can be useful when interpreting experimental data (l4 . 130T ] . 



IV. CONCLUSIONS 



The study presented here demonstrates the possibility of controlling the expansion of finite-size plasmas with 
appropriately shaped laser pulses. In fact, by irradiating a spherical target with sequential beams (double-pump 
technique) the phase-space dynamics of the ions can be manipulated, driving large-scale shock shells. These peculiar 
kinetic structures can lead to highly energetic ion-ion collisions and intraclustcr nuclear fusion reactions (in the case 
of clusters composed of deuterium or mixtures of deuterium and tritium). The features of the shock shells can be 
tailored by varying the intensity of the two laser pulses and the time delay between them, thus providing an efficient 
way to control the explosion of nanometer-sized plasmas. In the future, the availability of even more intense lasers 
will allow the control of the expansion of micrometer-sized, solid spheres, thus opening the way toward ultraprompt, 
ultralocalized sources of nuclear reactions, which can be useful for neutron production or nucleosynthesis. 

Since expansion control is achieved by acting on the electron dynamics, a detailed knowledge of the self-consistent 
dynamics of ions and electrons is required. The ergodic model described here provides a new, efficient tool, capable of 
describing the kinetics of spherical nanoplasma expansions with high accuracy. Although the ergodic model has been 
employed here in the paradigmatic case of initially Maxwellian electrons, it can also be used with arbitrary initial 
distributions of electron energy; this can be of interest in particular physical scenarios where accurate control over the 
expansion is crucial, such as in single-shot X-ray imaging of biological samples, in order to guarantee that no relevant 
sample damage takes place before the typical imaging time. 
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